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The measurability of the spin-orbit (S-L) couphng induced orbital effect is dependent on the 
orientation and magnitude of the orbital precession velocity, flo. This paper derives fto in the case 
that both spins in the binary system contribute to the spin-orbit (S-L) coupling, which is suitable 
for the most popular binary pulsars. Neutron star- White Dwarf star (NS-WD) binaries (as well 
as for NS-NS binaries). This paper shows that from two constraints, the conservation of the total 
angular momentum and the triangle formed by the orbital angular momentum, L, the sum the spin 
angular momenta of the two stars, S, and the total angular momentum, J, the orbital precession 
velocity, flo, along J is inevitable. Moreover, by the relation, S/L <C 1, which is satisfied for a 
general binary pulsar, a significant flo (in magnitude) is inevitable, 1.5 Post Newtonian order (PN). 
Which are similar to the case of one spin as discussed by many authors. However unlike the one spin 
case, the magnitude of the precession velocity of flo varies significantly due to the variation of the 
sum the spin angular momenta of the two stars, S, which can lead to significant secular variabilities 
in binary pulsars. 

PACS numbers: 04.80.-y, 04.80.Cc, 97.10.-q 



I. INTRODUCTION 

Barker & O'ConnellQ derived the first gravitational 
two-body equation with spins. In which the precession 
of the angular momentum vector, L, responsible for the 
S-L coupling is expressed as around a combined vector of 
the spin angular momenta of the pulsar, Si, and its com- 
panion star, S2 (S = Si + S2). Which is insignificant, 
2PN, typically 10^'^dcg per year, and therefore ignor- 
able in pulsar timing measurement. But notice that this 
ignorable orbital precession velocity, flo, is expressed as 
along a combined vector by Si and S2, which is not static 
to the line of sight (LOS). Therefore, the measurability 
of flo need to be further investigated. 

Since the precession of the pulsar spins (Si and S2) 
around L (and J) is significant, 1.5 PN (typically Ideg 
per year) by Barker & O'Connell two-body equation, and 
the precession of Si relative to the line of sight (LOS) has 
been measured through the structure parameters of pul- 
sar profile '5, "s*!, Therefore the small precession velocity 
of L with respect to a combined vector by Si and S2 
doesn't mean that it is also small relative to LOS. 

For binary pulsars the change of the orbital period 
due to the gravitational radiation is 2.5PN; whereas the 
geodetic precession of the pulsar is 1.5PN. Therefore, the 
influence of gravitational radiation on the motion of a bi- 
nary system can be ignored. And the total angular mo- 
mentum, J (J = L-t-Si-|-S2), can be treated as invariable 
both in magnitude and direction. Therefore after count- 
ing out the proper motion effect, J is at rest to LOS. So 
to make sense in observation, the precession velocity, i.e., 
flo, should be expressed as around J. 

In the study of observational effect of the orbital pre- 
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cession due to the coupling of the spin induced quadruple 
moment of the companion star with the orbital angu- 
lar momentum '3, 's', the precession of L is expressed 
as relative the total angular momentum vector, J. And 
the orbital velocity obtained can explain some secular 
variabilities measured. However this model might not 
suitable for neutron star-white dwarf (NS-WD) and neu- 
tron star-neutron star (NS-NS) binaries which have much 
smaller quadruple moment. 

On the other hand, Apostolatos et al0, KidderQ as 
well as Wex and Kopeikinj^ studied the precession of 
L as relative to J in discussing the spin-orbit (S-L) cou- 
pling in the NS-BH binary systems. Similar to the former 
quadruple-orbit coupling, this S-L coupling only consid- 
ers the spin (quadruple) of the companion star (BH), 
the spin of the pulsar is ignored. The difference is that 
the former Q-L coupling corresponds to 2 PN along and 
J, whereas the latter S-L coupling corresponds to 1.5 
PN along J. Therefore, the latter is easier to explain 
the measurements of NS-WD and NS-NS binaries, which 
have much smaller quadruple moments. 

The S-L coupling in the one spin case predicts a static 
precession of L, or flo is unchanged both in magnitude 
and direction. Therefore, the orbital effect induced in 
such case can be completely absorbed the orbital param- 
eters such as orbital period, Pb, and the precession of 
periastron, uj. However, this model cannot explain the 
significant secular variabilities measured in many binary 
pulsars. Therefore, it actually has not been used in pul- 
sar timing measurement. 

This paper studies the precession velocity of L around 
the static direction, J, in the a general case, two spins. It 
comes to the following conclusion: (a) as that of one spin 
case, the three vectors, L, J and S are still in the same 
plane at any instant; (b) as that of one spin case, the pre- 
cession velocity of L around J is significant (measurable 
for many binary pulsars); (c) unlike the one spin case, 
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the magnitude of the precession velocity of flo varies sig- 
nificantly due to the variation of the sum the spin an- 
gular momenta of the two stars, S, which can naturally 
explain the significant secular variabilities measured in 
binary pulsars 

This paper is arranged as follows. Section II introduces 
the orbital velocity of Barker and O'Connellllj and that 
of one spin equation by Kidder 8,], Wex and Kopeikin 6]. 
The measurability of them are also discussed and com- 
pared. Section III derived the orbital velocity as around 
J in the general (two spins) case. Section IV gives the 
derivatives of fJo, responsible for the significant secular 
variabilities. 



II. ORBITAL PRECESSION VELOCITIES BY 
BARKER & O'CONNELL AND OTHER 
AUTHORS 

According to the gravitational two-body equations 
the motion of the a binary system can be described by 
three vectors, the spin angular momentum of the pul- 
sar, Si, its companion star, S2, and the orbital angular 
momentum, L. The secular result for Si is 
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where ng , ni and W2 are unit vectors of the orbital angu- 
lar momentum, the spin angular momentum of the pul- 
sar and its companion star, respectively, mi, are the 
masses of pulsar and the companion star, respectively, e 
is the eccentricity of the orbit, is the semi-major axis. 
The first term corresponds to the geodetic (de Sitter) 
precession which represents the precession of Si around 
L, since L cx a^''^, the first term actually corresponds 

to a^^"^ (1.5 PN); and the second term represents the 
Lense-Thirring precession, Si around S2, which corre- 
sponds to oTj^ (2 PN). $72 can be obtained by exchanging 
the subscript 1 and 2 at the right side of Eq(^ . 

The precession of the orbital angular momentum, JIq, 
can be given as follows 
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where ilpAr is the relativistic periastron advance, 17^, is 
the geodetic precession caused by S-L coupling, and fig 
is the precession due to the Newtonian coupling of the 
orbital angular momentum vector to the quadruple mo- 
ment of the two bodies. For a general binary pulsar sys- 
tem, like NS-NS and NS-WD binaries, ilg is ignorable, 
2.5PN, SlpAT is significant, 1.5PN, but it has no contribu- 
tion to the S-L couphng, due to it is along L. Therefore, 
only J7i (2 PN) contributes to the S~L coupling. By 
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where a + 1 meant modulo 2 (2 + 1 = 1), EqQ means 
that L precesses around a vector combined by Si and S2 
with a very small velocity (2PN), and the precession of 
L due to geodetic precession is thus considered negligible 
to pulsar timing measurement. 

However since Si and S2 and thus the combined vector 
by Si and S2 precess rapidly (1.5PN) with respect to 
LOS, therefore from the small precession velocity L to the 
combined vector of Si and S2 (2PN) given by Eq@, we 
still cannot conclude that L also precesses with respect 
to LOS slowly. 

Actually in the study of the observational effect of the 
Q-L coupling, J7q is transformed to the direction of J. 
And the value of ilg along J is about L/S times larger 
than that of the previous one, thus ilg along J corre- 
sponds to 2PN instead of the previous 2.5PN. 

Similarly if orbital velocity corresponding to the S- 
L coupling is also expressed as along J, then it would 
also be L/S times the velocity $7^ given by Eq©, which 
would be 1.5PN instead of 2PN. Then it should be com- 
parable to ripAT in magnitude and measurable. 

Actually this work has been do by many authors || 
Q , in the special case of one spin. From the conversation 
of the total angular momentum, we have J = 0, which 
can be written, 



rjo X L = — rji X Si — X S2 



(4) 



Notice that as defined by Barker and O'ConnellU and 
Apostolatos et alQ, L = /zM^/^r^/^no, where M = m\ + 
7712, and = mxvnijM . 

In the special case that one spin angular momentum 
is ignorable, i.e., S2 — 0, as discussed by Kidder Wex 
and Kopeikin|^, (Apostolatos et al0 discussed the spe- 
cial case that Si • S2 = 0), S" is a constant in magnitude. 
And thus L and S precess about the fixed vector J at the 
same rate with a precession frequency approximately!^ 
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Note that Eq© uses G = c = 1 and r = ap(l - e^)^/^. 
Since L, J and S forms a triangle, therefore, the three 
vectors are in the same plane. Or more clearly, L and S 
are at the opposite sides of J at any instant (J is invari- 
able). 

Notice the velocity given by EqjS)) is relative to the 
static direction J (which makes sense in observation), 
and the magnitude of Eq© is 1.5PN, which means the 
precession velocity of L around J is much larger than 
that of L around a combined vector by Si and S2 given 
by EqQ. 

If Eq(|3Jl is applied to pulsar timing measurement, the 
static precession of L around J can be completely ab- 
sorbed by parameters like the orbital period, P;,, and the 
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advance of precession of periastron, uj. Therefore, the 
effects of Eq© might be unmeasurable, although it is 
significant. 

However, if both spins arc considered, then the pre- 
cession of Si and S2 will lead to a variable S (both in 
direction and magnitude), and therefore, through S-L 
coupling, the orbital velocity, fJo varies in magnitude. 
Which cannot be completely absorbed by parameters, 
like, Pb and Co. And it is this variation that automati- 
cally explains the secular variabilities measured in many 
binary pulsars 

Therefore, it is necessary to derive an orbital precession 
velocity that is relative to J and including two spins. 

III. ORBITAL PRECESSION VELOCITY IN 
GENERAL CASE 

A. The scenario of the motion of L, J and S 

This section derives orbital precession velocity as 
around J in the case that Si 7^ and S2 ^ 0, so that it 
can be applied to general NS-WD and NS-NS binaries. 
Eql^J can be rewritten as 

r2o X L = -r22 X S - (r2i - fla) X Si . (6) 

Notice that when ignoring the terms over 2PN, Jli and 
are along np. The first term at right side of Eq|(Bl) rep- 
resents a torque that is perpendicular to the plane deter- 
mined by L and S. And the second term at right side of 
EqlO also has a component in the same direction as the 

first one, which can be written as, (fii — ^2)S\ sinA^^^, 

in magnitude. Where s\ corresponds the component of 
Si in the plane determined by S and J, which satisfies 
s\ = Sicosrjssn and A|g^ is the misalignment angle 

between s\ and J. rjssi = VSi ~ Vs (jlSi and 775 are 
precession phases of Si and S respectively). 

Obviously S should precess at the velocity which is 
equal to that L precess around J, whereas, Eq® uses 
to replace this velocity, therefore, an extra term ((ili — 

f22)5'i' sin A|g^) is needed to make up the difference. The 
torques made by these two terms can be given, 

oil 

(Ti)fl = [f^asinALs + (f^i - ^^2)^ sinA^gJ^ , . (7) 

{ti)r is equivalent to a combined one, flsS sin A^s- The 
second term at right hand side of Eq® can produce a 
torque that is perpendicular to that of the first term. 
Which can be written, 

{T2)R^ini-n2)S^smXis^ , (8) 

where is the component of the Si, which is at the 
plane perpendicular to the plane determined by L and S 
instantaneously. Which produces a torque, (T2)fl, that 
is perpendicular to (ti)^, obviously in the one spin case, 



('''2)^ = 0. Actually Eq© can also be regarded as S pre- 
cesses around a vector that is perpendicular to the plane 
determined by L and S with velocity {t2)r/S instanta- 
neously. 

(''2)^ / (in the general case) indicates that the left 
hand side of Eql^J and Eq@ must make a correspond- 
ing torque to balance with (t2)_r, or the orbital precession 
velocity, ^Iq, much have a component that is perpendic- 
ular to the plane determined by L and S to balance with 
('''2)_R- Then it seems that the three vectors may not be 
in the same plane (or J is not in the plane determined 
by L and S) , due to {t2)r ^ 0. 

However this situation will never happen because the 
constraint, J = L -I- S, must be satisfied at any instant 
in the two spins case also. Which demands that J, L 
and S must be in the same plane at any instant. Si- 
multaneously, the constraint, J = 0, demands that J is 
unchanged both in magnitude and direction. Therefore, 
L and S have to be at the opposite side of J at any in- 
stant to satisfy to the two constraints at the same time. 
In other words, the plane determined by L and S rotates 
around a fixed vector, J. 

Then what is the role of the torque, {t2)r, in the mo- 
tion of the three vectors? It can cause a precession of S 
around a vector perpendicular to the plane determined 
by L, J and S (plane LJS) instantaneously, which can 
change the angle between L and S. Similarly L also pre- 
cesses around a direction perpendicular to the plane LJS 
instantaneously, which changes the angle between L and 
J. 

Thus {t2)r can only changes the misalignment angles 
of L and S with respect to J, it cannot influence the 
scenario that L and S precess rapidly at opposite side of 
J at any moment. In the one spin case {t2)r, = 0, the 
misalignment angles between S and L as well as L and J 
are unchanged. 

Therefore the difference between one spin and two 
spins is that the former corresponds to a simple preces- 
sion of L and S around J, but the latter corresponds to 
a nutation superimposed on an overall precession. While 
the common point is that at any instant, the three vec- 
tors are in the same plane, and moreover, L and S are at 
the opposite side of J. 

As discussed above from the two constraints, J = L-l-S 
and J = 0, we can conclude that the plane determined 
by L and S precesses around a fixed vector, J, rapidly 
(1.5PN). In other words, the precession velocity of L and 
S (J7o and fls) is along J. One can also choose vectors 
other than J as axis that L and S precess around. How- 
ever that would be like using a planet instead of the Sun 
as the center of a reference frame to explain the motion 
of the solar system. 

Actually the scenario of the motion of vectors in a bi- 
nary system have been discussed by authors, Smarr and 
BlandfordQ, Hamilton and Sarazin^^- From which one 
can easily deduce that the small velocity of L around a 
combined vector by Si and S2 doesn't influence that the 
three vectors, L, Si and S2 precess rapidly around J. In 
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other words, the precession velocity of L around J can 
be significant and therefore measurable. 

Since Barker and O'Connell's two-body equation sat- 
isfies the two constraints, J = L -I- S and J = 0. There- 
fore, Barker and O'Connell's two-body equation actually 
means that the orbital precession velocity, flo, is along 
J. However the misalignment angle between L and J is 
arbitrary in Barker and O'Connell's two-body equation, 
therefore, it is impossible to obtain Jlo along J without 
further constraint. 



B. The magnitude of fto 

Fortunately there is a constraint, S/L <C 1, that must 
be satisfied for a general binary pulsar. S/L 1 can 
impose strong constraint on X^j (misalignment angle be- 
tween J and S), and therefore, the magnitude (also di- 
rection) of Oq. 

By the triangle J = L + S and the relation S/L ^ 1, 
one can know that the misalignment angle between L and 
J must be very small, Xlj ^ 1- Thus no matter L par- 
ticipates in what ever motions, precession and nutation, 
L must always be very close to J. 

As discussed above, {ti)^ and {t2)r correspond to 
torques caused by components of Si and S2 that is in 
the plane LJS and perpendicular to it respectively. Then 
L must make a torque that equals (ti)^ in magnitude. 
Then L can react to (ti)^ in two ways: 

(a) L can precession around a vector combined by Si 
and S2 (which is the previous treatment of EqQ); 

(b) L precesses around J, with the opening angle of 
the precession cone A^j w S/{pL) {p = 1/ sinAjg). 

Mathematically, both (a) and (b) are allowed, (a) 
doesn't contradictory with (b), since that when L has 
certain angle with respect to the vector combined by Si 
and S2 (usually this angle is large), the angle between L 
and J can be arbitrary. Because the magnitude of L and 
S can be both comparable and very different {S/L <C 1) 
in Barker and O'Connell's equation. 

In (b), although the angle between L and J is very 
small, the angle between L and the vector combined by 
Si and S2 (or S) can still be arbitrary. Therefore, (b) 
doesn't contradictory with (a). The difference is that (b) 
uses one more constraint (which must be satisfied by a 
general binary pulsar) than (a), that is Xlj <C 1- Thus 
for binary pulsars, the orbital velocity obtained from (b) 
is more closer to the true one than that of (a). Moreover, 
the result from (a) has no observational means (not at 
rest to LOS), whereas the result from (b) make sense in 
observation (at rest to LOS). 

Therefore, for a general binary pulsar system, espe- 
cially NS-WD and NS-NS binaries, only (b) is considered. 
The left hand side of Eql^J and Eq® can be written as 

{ti)l = \^Qn J X Lno\ = floL sin Xlj ■ (9) 
in which Hq denotes the velocity of L around J. By Eq© 



and Eq®, we have the precession rate, 

^^0 = r ^^^l^ = pi}2 sin Xls + Pi^i -f^2)^sinA|5 . 
-L sm AL I o 

(10) 

Note pL sin Xlj ~ S. In reaction to the torque, {t2)r, of 
Eq®, L should precess around a vector, rij ( perpen- 
dicular to the plane LJS instantaneously), which corre- 
sponds to an opening angle of the precession cone tt/2 at 
any moment. 

(r2)L = ini^nj X Lno\ = n^.'^Lsin^ , (11) 

This precession is usually called nutation, superimposed 
on an overall precession of the system about the axis of 
total angular momentum. By Eq||SJ) and Ea ljlll) . we have 

^nu _ (^2)fl _ (»i - ^2)^1^ sin Afg^ 



Ea l|l()|l corresponds to a rapid precession of L around J 
(1.5PN); whereas Ea H12|) corresponds to a slow precession 
(nutation) of L around a vector that is perpendicular to 
the plane LJS at any instant (2PN). In other words, the 
velocity of Ea (|12l) derived from the torque, {t2)r, can 
only change the misalignment angles, Xlj and Xls, it 
cannot influence the direction of SIq. 

In summary, from the constraints, L + S = J, and 
J = 0, the precession velocity flo along J is inevitable. 
Further by the constraint, S/L <^ 1, which is correct for 
a general binary pulsar, a significant flo (comparable to 
Lu the precession of periastron) is inevitable. 

Therefore, the one spin case and the general case are 
similar in serval aspects. The orbital precession veloci- 
ties of a binary pulsar system are both significant, 1.5PN 
and both around J. Moreover, the three vectors, L, J and 
S are in the same plane at any instance in both cases. 
The one spin case corresponds to a zero nutation veloc- 
ity, while the two spins case corresponds to an ignorable 
nutation velocity of the orbit plane, 2PN. 

However, there is significant difference between one 
spin and two spins cases, the former corresponds to a 
constant S, which means a static precession of S and L 
around J, = const; whereas, the latter corresponds 
to a variable S, and therefore, a variable fig, which is 
responsible for the significant secular variabilities mea- 
sured in many binary pulsars. On the other hand, the 
significant secular variabilities measured in many binary 
pulsar seems to support the two spins case|^. 

IV. DERIVATIVES OF 

By Ea l(Tn|l . if S is unchanged then L will precession 
with a static velocity, flo, around J. However, since Si 
and S2 precess with different velocities, ili and ll2, re- 
spectively. Thus S varies both in magnitude and direc- 
tion (Si, S2 and S form a triangle), and in turn Xlj 
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and Xls elIso vary with time (S, L and J form a triangle, 
Fig.l). 

The motion of the vectors Si, S2 and S can be studied 
in the coordinate system of the total angular momentum, 
in which the z-axis directs to J, and the x- and y-axes are 
in the invariance plane. S can be represented by and 
, the components parallel and vertical to the z-axis, 
respectively: 



5 = (^^ + 5^)1/2 
and can be expressed as 

~ Si cos Xjsi + S2 cos A 



(13) 



JS2 ; 



= [{Srr + (S^r - 2S^S^ cos ,7s, sj 



11/2 



(14) 



SY = Si smX jsi and S2 = iS'2sinAj52 represent com- 
ponents of Si and S2 that are vertical to the J. SY , SY 
and form a triangle. the misalignment angle 

between SY and 5*^ can be written 



VS1S2 = (^1 - ^^2)^ + 00 



(15) 



where 



Xi = -noiXita,Il1]sSi ~ (f^l2^2^3 



X2 = -^01 sec^ rissi 



X3^Ya + Y& 
where Y = , a 



a sin Xjs ' 
2, 



Y = ni2{YcotTjs,S2-Y'a) 



[sin XjsXi 



(asinAjs)2 



-F(2 sin A js +^4) 



cos^ Xjs 
sinAjs 



sin Xjs 



[X3X4 + X4X3 



X^flu 



2VI12X3X2 



sin A 



JS 



Note that S, Xls-, S\ {S\ — Sicosrjssi), and A^^^ in 
Eg 1)10(1 are function of time, which leads to the change of 
the precession rate of orbit, 



X4 — —^112X3X4(2 



X4 



sin A 



LS 



(16) 



= Xls, cot A^s^ - {X[sf csc^ X[s^ 



where 

f2 = ^2^112X3X4 — ^12X1 (rioi^2 + ^V2^?>) + ^12^^1X5 , 

(17) 

where r2i2 — ^\ — ^2-, f^oi = f^i — i^o, 



Xi 



s 



sin XlSi : ^2 = tan?7ssi , X3 



SYSY sin ^,1,2 
a sin Xjs 



where A^^^ = tan ^ (tan A js^ cos 7755 J , A|g^ 



A 



z(l+z-')-2z 



LSi - — (T+i'p — ■ which z = tanAjs^ cosrjssi , 
z = -JloitanAjSi sinryssi , 2 = -Sloi tan Ajg^ cos ryss^ 

By Ea (|18|l . the third order derivative of fio are given, 



cos^ Xls ^ 



sin A 



cot A r c A 



X'3X4^-\ 



LS 



LSi'^LSi ' P 



sin A 



JS 



witha = sinA,5 + f^. 

Note that f2i and f22 are unchanged, and Als„ are 
unchanged (since they decay much slower than the orbital 
decay.?]. By Ea (|16|l . the second order derivative of JIq 
are given. 



^Iq^ pVL + 2pVl + pVt , 



(18) 



where 



57 — Vl2^i2{X3X4 + X3X4) — VIi2Xi{VIqiX2 + $112^3) 



where 



.... 

— — - — 5725712(2X3X4 + X3X4 + X3X4) — 



25712^1 (5701^2 + 5712X3) — 57i2Xi(rioiX2 + 5712X3 



^^12X1(5791X2 -l- 5712X3) -t- rii2(XiX5 -I-X1X5 -|- 2X1X5) , 

(21) 

Xi = -57oiXi tan?75Si - ^liXi sec^ 77551 



-5712X1X3 — 5712X3X1 



5^12X1(5701X2 + 5712X3) + 57i2(XiX5+XiX5) , (19) 



X2 = 257o^ tanr/ssi sec^ 7755^ 



X3 = Ya + aY + 2Y& ^ = -(2 cos Xjs + cos^ Xjs csc^ Xjs)Xjs 

X4 = -ni2{2 + COtXLs){X3X4 + X4X3)- = ^12X3 COS Ajs CSC A JS 

ni2X3X4{X4 CSC XlS - COS^ XlS CSC^ XlS^12X3X4) y _9^\ll ^2 x|| 3 



^5 = 2(A:^sJ^cosA:^^^ csc^A:^^^- 



where 



F = Ql2Fc0t7?s,s.-^^?2^CSc2 775,5, -4ni2yF(7-2ni2F2a .|| ..|| 2,11 ,J| 

'^■^LSi'^LSi ^LSi "I ^^3~ cot A^5j 

a = 2cTCT^ (a sin Ajs)^ + cr^ [sin AjsX4 + X4A7S COS Ajs d^A" (i+«2)^_2222-2i3-4(i+22)2iA" 

where = ^ ^^^p 

^ = Ogi tan AjSi sin r]ssi ■ 

+12i2 (2 A js COS Aj5 + X4) cos^ Aj5 CSC AjsXa 

+Or2(2sin A,5 + X4)iiX3 + ^Xa)] Acknowledgments 
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FIG. 1: Angles and orientation conventions relating vectors, 
S, L and J to the coordinate system. S, L and J form a 
triangle, and determine a plane LJS. At any instant S and L 
must be at the opposite side of the fixed vector, J (invariable 
both in direction and magnitude). and Si are compo- 
nents of Si projected in the plane LJS and perpendicular to 
it respectively. {ti)r and {t2)r are torques corresponding to 
the precession and nutation of L (or S) in a binary system 
respectively. 



